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ABSTRACT 

As a formation route for objects such as giant planets and low-mass stars in protostellar discs 
(as well as stars in AGN discs), theories of self-gravitating disc fragmentation need to be 
able to predict the initial masses of fragments. We describe a means by which the local Jeans 
mass inside the spiral structure of a self-gravitating disc can be estimated. If such a self- 
gravitating disc satisfies the criteria for disc fragmentation, this estimate provides a lower 
limit for the initial mass of any fragments formed. We apply this approach to a series of self- 
gravitating protostellar disc models, to map out the typical masses of fragments produced by 
this formation mode. We find a minimum fragment mass of around 3 Jupiter masses, which is 
insensitive to the stellar mass, and that - within the parameter space surveyed - fragments with 
masses between 10 and 20 Jupiter masses are the most common. We also describe how the 
Jeans mass allows us to derive a more general criterion for disc fragmentation, which accounts 
for the processes of viscous heating, radiative cooling, accretion and the disc's thermal history. 
We demonstrate how such a criterion can be determined, and show that in limiting cases 
it recovers several fragmentation criteria that have been posited in the past, including the 
minimum cooling time/maximum stress criterion. 
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1 INTRODUCTION 

During the process of low mass star formation, an excess of angu- 
lar momentum in the parent molecular cloud will generally result 
in the formation of a circumstellar disc. In the early phases of this 
disc's existence, it is expected to be relatively massive, cool, and 
extremely weakly ionised. Therefore, disc self-gravity is expected 
to play a dominant role in protostellar and protoplanetary di sc evo- 
lution JLin & Pringlell987l : lLaughlin & Bodenheimeill994h . in the 
same manner it is expected to play a role in the evolution of Ac- 
tive^ Galactic Nucle i (AGN), which also possess accretion discs 
JShlosman & Begelmanll 19891 : iGoodmar 12003). A self -gravitating 
disc becomes gravitationally unstable if JToomrdl 19641) : 



Q 



ttGE 



(1) 



where Cs is the local sound speed, S is the surface density and 
K is the epicyclic frequency (for Keplerian discs, this is equal to 
the angular frequency Q,). This result applies primarily to axisym- 
metric perturbations, and is usually extended to non-axisymmetric 
perturbations by changing the critical Q value to Q < 1.5 — 1.7, 
a result established empirically by numerical simulations (see 
IPurisen et al.ll2007l for a review). 
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The onset of gravitational instability (GI) produces non- 
axisymmetric, spiral structures. It has been shown that self- 
gravitating discs will typically settle into what is known as a 
marginally stable state, where the stress heating produced by the 
spir al waves is balanc ed (at least approximately) by radiative cool- 
ing jPaczvnskjl 119781) . In the marginally stable state, the spiral 
waves become quasi-steady structures that, while individually tran- 
sient, as a whole remain for long timescales, mediating the trans- 
port of angular momentum outward (and the transport of mass in- 
ward). 

The evolution of self-gravitating discs can therefore be neatly 
parametrised using the phenomenology of viscosity, much in 
the same way that the action of magneto-rotational instability 
(MRI) has been parametrised in less massive, mor e ionised discs 
JBlaes & Balbu3(l994l : Isalbus & PapaloizoJI 19991) . This pseudo- 
viscous parametr isation, and its encapsula tion in the dimension- 
less a parameter JShakura & Sunyaev 1973), has been a useful tool 
in si mplifying the computation of the structure of self-gr avitating 
discs teafiko\^l2005l : IClarkel2069l : lRice & Armitagell2009t), as well 
as co mputing the ev olution of s aid discs fLvnden-B ell & Pringld 
1 1974 Pringle 1981: L in & Pring le 1990; Rice et al. 20i^^ 

While useful, this pseudo-viscous approximation is precisely 
that - an approximation, which is accurate only for a limited re- 
gion of the available disc parameter space. In particular, for the 
approximation to be valid, the angular momentum transport that 
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occurs must be local. Gravity is an inherently non-local force, 
and if the spiral structures in the disc can exert sufficient gravi- 
tational field stresses on fluid elements at large separations, then 
the pseudo-viscous approximation begins to fail, and the disc's 
angular momentum transp ort is governed by global properties 
JBalbus & Papaloizoulll999i) . The breakdown of the local approx- 
imation has been shown to occur when the disc-to-star mass ra- 
tio, q, is gr eater than 0.5, and/or the disc scale h eight satisfies 
H/r > 0.1 JLodato & Ricel2004l : lForgan et alJlOlU) . Discs which 
satisfy these conditions typically possess strong m — 2 spiral 
structures, which exert significant torques from large distances, and 
exhibit highly variable temperature structures, and consequently 
highly variable Q. We will see that the issue of whether disc an- 
gular momentum transport is local or global will come to be of 
great importance in this paper. 

If the transport is indeed local, the balance between the magni- 
tude of disc stresses and the magnitude of radiative cooling has also 
been useful in describ ing the fragmentation of local self-gravitating 
discs. iGammid J200 11) demonstrated that a second criterion, along- 
side the Q-criterion above, must be satisfied if a self-gravitating 
disc is to fragment. This criterion is related to the local cooling 
time, tcooi, and is usually expressed using the dimensionless cool- 
ing time parameter /3c = tcooi^, i.e. 



Pc ^ Pcrit 5 



(2) 



where the exact value of Peru is typically established by numerical 
simula tions. Using the 2D shearing sheet approximation, lOammiel 
( 120011) showed that when the ratio of specifi c heats in two dimen- 
sions 72C — 2, Pcrit = S. lRice et al.l ( I2OO3 ) extended this analysis 
to 3D smoothed particle hydrodynamics (SPH) simulations (where 
Pc was fixed at some value), and explored the dependence of Peru 
on 7. They were able to confirm Gammie's assertion that, in the 
case where the disc has reached thermal equilibrium (i.e. the radia- 
tive cooling and stress heating due to gravitational instability are in 
balance), the fragmentation criterion can be recast in terms of the 
a parameter: 



(3) 
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where the critical value of a = acru ~ 0.06 for all values of 7. 
The fragmentation criterion can now be interpreted as a maximum 
GI stress that the disc can support in thermal equilibrium without 
fragmentation. 

The inner regions of discs are typically too hot for 
Q ~ 1, and are therefore not gravitationally unstable and 
incapable of fragmentation. Also, if the primary contribution 
to a is from GI stress, it will typical ly decrease w i th in- 
creasing proximity to the centra l star JArmitage et al.l uOOll : 
iRice & Armitagd l2009l : iRice et al.l l2010h . These facts immedi- 
ately suggest a minimum radius at which fragmentation can 
occur. Many authors, using a variety of techniques, including 
semi-analytic methods, grid-based and particle based simula- 
tions, find that for self-gravitating protostellar discs the minimum 
is ap p roximately 40 au, using typical disc parameters JRafikovl 
20051: iMatzner & Levirj l200a IWhitworth & StamatellosI I2OO6I : 



Bolev et alj 12004 IStamatellos & WhitworthI I2OO8I : iForgan et all 
20091 : IClarkell2009l) . 



This also constrains the maximum available mass to form 
a fragment (as typically equilibrium surface density profiles de- 
crease with increasing distance). Despite such constraints, the con- 
sensus view is that the objects formed occupy the higher mass 
end of the planetary regime (and the low mass end of the stellar 



regime). iKratter et al.l 1 1201 Of) estimate the available mass within 
the most unstable wavelength to fragmentation (the disc scale 
height H) using one-dimensional disc models to obtain fragment 
masses, and then model the fragments' growth embedded in the 
disc to constrain expected planet masses. Studies o f the behaviour 
of Pc (or more strictly tc poi) in similar disc models |Rafikovll2005l : 
iNero & Biorkmanll2009l) tend to agree with lKratter et al.l ( I2OI0I) in 
that they calculate a minimum fragment mass of around 5Mjup, 
where Mjup is the mass of Jupiter 

While the Q and Pc criteria provide an elegant means by 
which to predict disc fragmentation, there are complicating factors 
that must also be considered. Perhaps the most important is the fact 
that (for the most part) these numerically simulated discs are iso- 
lated. We expect that self-gravitating discs - as they form at an early 
phase in protostellar evolution - will be surrounded by a substantial 
envelope from which the disc can accrete. The potential for mass- 
loading as a result of accretion can alter the local opt ical depth as 
well as the local stresses felt by the disc. [Kratteret_alJ (2010) show 
in grid-based numerical simulations that systems accreting rapidly 
can sustain a ^ 1 w ithout fragmenting. This is also confirmed by 
iHarsono et al.l ( 120111) using SPH simulations, where they ascribe 
the extra Reynolds stresses to velocity shear at the disc surface. 
Both show that the disc infall rate and disc-to-star accretion rate 
are typically proportional to each other, and that non-local angular 
momentum transport is significant. These factors alone show the 
limited validity of isolated, constant /3c simulations in estimating 
disc fragmentation criteria. 

Further to this, the imposition of a fixed Pc is somewhat un- 
physical. As the disc cools through varying opacity regimes, the 
opacity is expected to ch a nge significantly (s ee e.g. iForgan et al.l 
'200^: ICossins et al.l2009l) . IClarke et al] ( l2007h carry out SPH sim- 
ulations where Pc is initially large, and gradually decreased towards 
the critical value. They find that the eventual Peru can be up to a 
factor of 2 lower This reduction of PcrU wil l incre ase Ocru by the 
same factor - while the results of iRice et al.l ( 120051) provide a use- 
ful watermark, the exact va lue of OcrU is not precisely determined. 
Indeed, a related result of ICossins et al.l ( 12010) shows that Peru 
also depends on the opacity regime the gas is in, producing isolated 
discs which require aeru = 0.1 for fragmentation. 

Taking these factors into consideration, it becomes clear that 
the current use of Pc or a is not sufficiently general as a fragmenta- 
tion criterion. While it is suitable for self-gravitating discs of con- 
stant Pc that do not accrete, such discs are not typically found in 
astrophysical situations. Ideally, we would like to move towards a 
situation where the fragmentation criterion is sufficiently rich in its 
content that it can describe the full complexity of the fragmentation 
process in realistic self-gravitating discs. 

This paper introduces two useful analytical formalisms: the 
first describes the local Jeans mass in a marginally stable, self- 
gravitating disc (which is described in section |2j; this provides a 
means of predicting the masses of objects formed as a result of disc 
fragmentation. In section|3]we demonstrate that this expression for 
the Jeans mass can be used to develop a new timescale criterion, 
which folds in the influences of stress heating, radiative cooling, 
accretion and thermal history in general. Having developed these 
concepts, we construct simple one-dimensional models to investi- 
gate the typical fragment masses produced by self-gravitating discs 
(described in section|4), the results and discussion of which can be 
seen in section[5] We summarise the work in section|6] 
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2 THE JEANS MASS INSIDE A SPIRAL WAVE 
PERTURBATION 

While it is perhaps more appropriate to look at the Toomre mass in 
this rotating configuration, the Jeans and Toomre masses are rel ated 
by fac tors of order unity when the disc is marginally stable ( Nel sonI 
|200g), so it is safe to use the more conceptually simple Jeans mass. 
As fragments will typically form thanks to the perturbing influence 
of a spiral density wave, it is sensible to look at the Jeans mass 
inside the wave itself. 

Let us assume that we have a marginally stable disc (Q ~ 1), 
with some local surface density E, producing spiral wave perturba- 
tions with fractional amplitude ■^. Essentially, we want a pertur- 
bation to contain a mass M greater than the local Jeans Mass, Mj. 
If we assume a spherical Jeans Mass: 



M,J = -TT 

3 



Gp, 



■pert 



3/2 



Ppert — ^ TT 



5/2 



G^/^p: 



1/2 ' 



(4) 



where Ppert is the density of the perturbation, and Cs is the local 
sound speed. In the midplane. 



Ppert — ^pert/^-^^ 

where by definition 



^pert — ^ 

Substituting into Mj gives 



-f) 



IH^/^ 



_ 4:V2 5/2 

■' ^^~'' G3/221/2 (1 + A^) 

We now eliminate E using Q: 



Nl/2 



Csfi 



(GE)-^^^ 

This gives 

4^7r^Qi/2 



1/2 



Mj^ 



5/2^1/2 



3G 



ny^i + ^) 



Using H = Cs/Q,, we can eliminate Q to obtain: 
4^271^ Qi/2c2^ 



Mj 



3G 



1 + 



(5) 



(6) 



(7) 



(8) 



(9) 



(10) 



This is a general description of the Jeans Mass inside the spiral 
wave, achieved using the following assumptions only: 

(i) The disc is marginally stable (Q ~ 1) 

(ii) p = T,/2H at the midplane 

(iii) H = Cs/n 

(iv) The Jeans mass is spherical 

Using the non self-gravitating prescription for H might seem haz- 
ardous, however for marginally stable discs 



J^ns 






Hs 



ttGE' 



(11) 



(this can be confirmed by substituting for Q ^ 1). Assumption 2 is 
reasonable for cases where the thin disc approximation applies (as 
is the case for assumption 3). We can calculate the fractional am- 
plitude AE/E using the empirical results of iCossins et a l. (2009): 

AE _ 1 



(12) 



This now gives 



Mj = 



aV^tt'^ Q^ 



3G (1 + 1/^,) ■ 



(13) 



This quantity is now completely calculable using azimuthally aver- 
aged variables. The Jeans mass depends on the local sound speed, 
the angular velocity and the thermal physics of the gas (which is 
sensitive to the surface density and temperature structure). We must 
however note that using this empirical prescription for the frac- 
tional amplitude requires that the angular momentum transport in 
the disc be local. 



3 A GENERAL TIMESCALE CRITERION FOR DISC 
FRAGMENTATION 

Can the Jeans mass provide a timescale criterion for disc fragmen- 
tation (much as the cooling time has done in the past)? Let us now 
define three timescales, normalised to the orbital period: 

Mj 



-n, 



E 



(14) 
(15) 
(16) 



We use /3 and T to distinguish between variables of differing be- 
haviour. Pc is positive definite - this is to keep it in line with the 
conventional definition as it is currently used (and also to pre- 
vent AE/E taking imaginary values). Fs measures the competi- 
tion between disc accretion and stellar accretion. Vj measures the 
timescale on which the local Jeans Mass changes - in Vj orbital 
periods, we can expect the Jeans mass to either double or decrease 
to zero, depending on the sign of F,j. E and Mj can be either posi- 
tive or negative - for fragmentation to be favourable, the local Jeans 
mass must be decreasing. A small, negative V j therefore represents 
the most likely circumstances for disc fragmentation. The critical 
value of V J is less clear - any disc that can maintain a negative T j 
will proceed towards fragmentation. Rather than presenting us with 
two discrete regimes, fragmenting and non-fragmenting, we see a 
spectrum of potential outcomes, some fragmenting rapidly; some 
fragmenting on much longer timescales (possibly longer than the 
lifetime of the disc, and therefore effectively non-fragmenting); and 
others moving away from potential fragmentation (either slowly or 
rapidly). We suggest that — 10 < Fj < gives highly favourable 
conditions for prompt fragmentation, but the lower limit is by no 
means fixed, and will require empirical confirmation. 

We wish to now derive F j: substituting for the self- gravitating 
scale height, we can rewrite the Jeans mass in terms of E, Cs and 



Mj = 



4^271^ 



3G2 E(l + l/^/:fl; 



To derive Mj, we can use the chain rule: 

dMj 



Mj 



dMj . , dMj- 
-Cs H — rr^^E ■ 



9cs 



aE 



dp. 



-Pc 



(17) 



(18) 



Instead of directly calculating Cs, it is easier to describe the rate of 
change of specific internal energy u, which we can relate back to 



7(7- 1)m. 



(19) 
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We will assume that the sound speed varies due to two effects only 
- firstly, radiative cooling, given by 

u 

and secondly, heating due to viscous dissipation, given by 



(20) 



Uheat = 9/4q7(7 - 

Using the chain rule on Cs and rearranging gives: 



dcs 



i/2cs 9/407(7 - i)n 



1 



(21) 



(22) 



CLU y T-cool 

We now use equation (TTJl to evaluate equation (|T8), where we sub 
sequently substitute back for Mj: 



Mj 



4 c. f 9a^i-f-l)n 

^^^T i 4 



1 

T'cool 



Mj 



+ 



M 



J 



40^ (1 + l/^/^) ' 
Equivalently, 

1 



(23) 



M. 



+- 



MjQ I 2 



13c 



+ 



907(7 — 1) 



1 



This quickly gives the Jeans timescale Fj as 

/?c ^ 4 ) ^ r^ 



(24) 



r,= (2(-l + 5^3ll^ 



/3c 



Ap^J^n (1 + i/VP^) 



(25) 



Remember that we require the magnitude of Tj to be small - 
making it more negative means that the Jeans mass will decrease 
on a slower timescale, making fragmentation less likely. Also, 
small positive values of Tj will rapidly increase the Jeans mass, 
suppressing fragmentation. 

To recap, the following assumptions have been made: 

(i) The disc is marginally stable (Q ~ 1) 
(ii) p — T./2H at the midplane 
(iii) H = c,/n = cI/ttGT, 

(iv) AE/E = 1/v^ (which implicitly assumes locally trans- 
porting discs) 

(v) The Jeans mass is spherical 



3.1 The Limiting Values of Tj - Recovering Older Criteria 

3.1.1 The "Standard" Cooling Time Criterion 

If we assume E = 0, then -^ >■ 0. If we also assume Be = 0, we 



r. 



Pc 



+ 9/407(7 - 



(26) 



and the fragmentation is governed by the cooling rate (or equiv- 
alently the disc's GI stresses). The minus sign appears due to our 
definition of /3c : this is sensible considering that cooling will reduce 
the local Jeans Mass. If we substitute for a in the case of thermal 
equilibrium 



97(7 - l)/3c 



(27) 



then the cooling and heating are in detailed balance and Tj — )■ 00, 
i.e. AIj = const.. This is essentially a re-statement of the assumed 
thermodynamic equilibrium, and is almost axiomatic. If the grav- 
itational torque saturates at some maximum possible value of a, 
then continually increasing the cooling rate will eventually force 
r.7 to become small and negative, rapidly decreasing the local Jeans 
mass until fragmentation becomes possible even for relatively low 
density perturbations. We must once more note that local, pseudo- 
viscous transport is assumed. If the disc angular momentum trans- 
port is global, we must replace the heating term with a more com- 
plicated expression reflecting non-local effects. Indeed, even discs 
with local angular momentum transport will be subject to non-local 
heating such as from stellar irradiation - this extra heating term 
should also be considered if modelling/simulation permits. 



3.1.2 The Effect of Thermal History 

If we allow Pc to be non-zero, then we can compare two limiting 
cases. If the change in cooling time is low, then we are in the regime 
described in the previous section. If however the change in cooling 
time is rapid: 



Pc » Pc^ 



4n/3g''' (1 + l/v%) 

Pc 



(28) 



with r ,j becoming small. This is in accordance with lClarke et al.l 
( 120071) , who show that discs with low Pc are generally more stable 
to fragmentation than those with high (negative) Pc. But what is 
responsible for the disc's fragm entation? We can investigate this 
by applving lClarke et alj ( l2007n 's prescription for Pc{t): 



Pc{t) = PdO) 1 



t < T 



where T defines the magnitude of Pc'. 

PM 



Pc 



T 



(29) 



(30) 



Substituting into equation ([25}, this gives 



r.7 = 



/?c(o)(i-i/r) 
1 



+ 



907(7 — 1) 



4™(V/3c(0)(l-t/T) + l) 



(31) 



Let us consider two cases: if there is no maximum stress beyond 
which fragmentation occurs, then the first two terms will always be 
able to balance, Pc will eventually tend towards zero by definition, 
and 



-4™ 



(32) 



If Pc is large, T will be small and hence Fj will be small and 
negative, promoting fragmentation. If there is a maximum stress 
which the disc exceeds, thermal equilibrium is no longer possible, 
and Fj will continue to decrease towards small values. If we fix 
the viscous heating term at its maximum, we can recover the rate 
of change of Fj against time in this scenario: 



dFj _ p^jPcjO) 
dt TO, 



_1_ 



8^fi^/A^(^/A^+l) 



(33) 
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The derivative is positive, and its magnitude depends primarily on 
tiie rate of chiange of j3c and on T j itself. For discs that are rapidly 
changing their cooling time, T is small, and hence Tj will increase 
rapidly towards fragmentation in only a few orbits. For discs chang- 
ing their cooling time on longer scales, then T and T ,j are both ini- 
tially large, and therefore the derivative is also initially large. How- 
ever, as the derivative decreases rapidly with decreasing F,/, large 
T discs will proceed much more slowly towards fragmentation, re- 
quiring many more orbits. In fact, the number of orbits required for 
fragmentation in this case will depend linearly on T - increasing 
T by a factor of 100 will require 100 times as many orbits to wit- 
ness fragmentation. If T is sufficiently large, then secular evolution 
of the disc (e.g. significant stellar accretion) can act to reduce the 
likelihood of fragmentation by removing mass from the system. 

A similar result applies to simulat ions where 13c is re duced and 
then held at some value Puoid (such as lClarke et al.l2007h . This ac- 
tion will fix V,j at some value given by the mismatch between the 
first and second terms in equation M5\ . Low values of Phoid will 
increase the mismatch, decreasing Fj and reducing the timescale 
for fragmentation to occur. Given this and equation ( |33t, we can 
now m ake some general statements about the work of lClarke et al.l 
( 120071) . They see some discs fragmenting at very low /3c because 
once favourable fragmentation conditions arise, there is a finite 
timescale for fragmentation to occur, defined by V ,j. If T is suffi- 
ciently small (to the point of being unphysical), then j3c can become 
extremely small before T j reaches values necessary for fragmen- 
tation. Also, we show that as C larke et al. (2007) note themselves, 
if they were able to run their high T simulations with /3c below the 
critical value for a sufficiently long time, they would also proba- 
bly fragment. The length of simulation time required to produce a 
fragment will be calculable from the evolution of F,j. While high 
T simulations may be more stable to fragmentation in the short- 
term, their long-term stability is not likely to be much better than 
simulations where the cooling time is rapidly changing. 

As a final aside on the subject, we must be careful about the 
sign of /3c. Fragmentation is more likely if T j is negative, and 
hence /3c should be negative also. If the change in cooling time 
is rapid and positive, then /Jc will become too large, and the Jeans 
mass will either increase or become steady at a higher value (which 
can be seen from the /Jc dependence of equation MIV ). Large fluc- 
tuations in /3c are therefore not a guarantee of fragmentation - we 
require the cooling rate to increase rapidly for fragmentation to be 
more favourable. This is akin to the fragmentation of shocked flows 
producing star formation in molecular clouds (e.g. iHeitsch et al.l 
[2008; Bonnell & Rice 2008), and may explain the formation of ob- 
jects during protos t ellar encounters w i th sufficiently exte nded discs 
dThies et al.l l2005l : IShen et alj I2OIOI: iThies et alj l2010h . Equally, 
this accounts f or fragmentation suppression during enco unters with 
compact discs ( iLodato et al.l20 07: Forgan & Ricel200 9). where the 
subsequent increase in optical depth ensures /Jc is large and posi- 
tive. 



material does not merely increase the disc mass - it affects the local 
thermodynamics and the angular momentum distribution, both of 
which can push the disc away from Q ^ 1. We should therefore 
also demand 



dQ_dQ_dc^ dQ_dn dQ ^ 
dE ^ dcs dE "*" 90 dE "^ 9E ~ 



(35) 



if the disc is to remain amenable to fragmentation. Substituting 
dCs/dT, — Cs/T, (and similarly for fi), we obtain 



dQ ^ Qs^,Q^_Q <(. 

dE^c, Ei^E S~ 

Upon rearrangement, the following condition now appears: 



(36) 



Cs n 



<o 



At the critical point (where -^ 
merely 

E ^ c, n 

E ~ cs n' 



(37) 
0), this can be simplified to 

(38) 



If accretion results in greater deposition of angular momentum or in 
local heating, then the disc moves away from fragmentation (this is 
intuitively obvious from the definition of Q). Also, it further under- 
lines the stabilising influence of sus tained heating, either by accre- 
tion shocks or by stellar irradiat ion ( lRafikovll2007l : ICai et al]|2008l; 
IStamatellos & Whitworthll2008h . 



3.1.4 Summary 

Marginally stable self-gravitating discs are susceptible to fragmen- 
tation if: 

(i) The cooling rate is sufficiently high (or equivalently the disc 
gravitational stresses are sufficiently high), or 

(ii) The cooling rate increases sufficiently rapidly (or the disc 
gravitational stresses increase rapidly), or 

(iii) The disc's envelope accretion rate is sufficiently high (while 
maintaining Q ~ 1). 



3.1.5 Predicting Fragmen t Masses 

We should be able to identify an equilibrium Jeans Mass by de- 
manding Mj = 0, or 



1 



+ 9/407(7 - 1) 



1 



+ 



Pc 



0(39) 



/3c ' "'' "' '' Te ■ 4/3c'^2n (1 + 1/V^^ 

This is a first order non-linear ODE for /3c (assuming a constant 
Fe). In the constant /3c formalisms of the past, this becomes 



3.1.3 Accretion as a Shortcut to Fragmentation 

Let us assume for the moment that heating and cooling are in exact 
balance (with a large, constant value of /3c), such that the heating 
and cooling terms in equation|25]cancel and /3c — 0, and hence we 
can approximate 



Tj = -Fi 



(34) 



i.e. glibly speaking, discs which have a higher accretion rate are 
more favourable to fragmentation. Of course, accretion of envelope 



/3c = 9/407(7 - 1) - Fe 



(40) 



These formalisms are not predictors of the fragment mass per se, 
but they show the secular Jeans mass, and to some extent the secular 
fragment mass if the disc is left to its own devices and satisfies the 
correct criteria. Also, in practice the cooling rate will be a function 
of deposition rate (in the sense that this maintains a local optical 
depth), the angular velocity will be a function of mass ratio and the 
cooling rate will be sensitive to changes in temperature and surface 
density. This differential equation is only soluble numerically, and 
probably only with numerical simulations. 
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In the case of simple analytic models, it is more expedient to 
impose the fragmentation criteria, and calculate M,j at the frag- 
mentation boundary. This method will now be applied to such disc 
models in the following sections. 



4 FRAGMENT MASSES IN SIMPLE DISC MODELS 
WITH LOCAL ANGULAR MOMENTUM TRANSPORT 

To investigate the dependence of the Jeans mass on disc parameters, 
we construct a model of a self-gravitating disc assuming a steady 
state accretion rate M and an outer di sc radius rput . We take an ap- 
pro ach essentially identical to that of lLevinI J2007l) . IcIaA^ ( |2009|) 
and lRice & Armitagel ( 120090 . For the sake of simplicity, we assume 
E = 0, /3c is small, and impose — 5 < T j < as the criterio n for 
fragmentation. Cautioned by the results of lCossins et al.l ( l2010l) , we 
impose a saturating value for a = asat = 0-1. 

The model is constructed as follows: we assume the disc is 
Keplerian and gravitationally unstable at all radii, i.e. 

for any value of r. The steady state accretion rate is 

37rQC„E 



M 



Q. 



(41) 



(42) 



which is assumed to be constant over all radii. We will also assume 
that a is determined by the constraint of local thermal equilibrium, 
i.e. 

4 



(43) 



97(7-l)/3c' 

where /3c = {u/ii) f2 is a function of volume density p and temper- 
ature. This assumes of course that the disc is local - we will leave 
a discussion of the validity of this a parametrisation for a later sec- 
tion. The cooling function is 



Tn4 

OSBl 

r+1/7 



(44) 



where asB is the Stefan-Boltzmann constant, and the optical depth 
T = 'Sk{p, Cs), with K the opacity. Our aim is to calculate the 
three unknowns a, Cs and E at any radius. Equations J41b . J42b and 
( 143 1 place three constraints, and the system is therefore soluble. We 
employ the equation of state of Stamatellos et al. (2007) (see also 
iForgan et al.ll2009l : ICossins et al.ll2009l : iRice & Armi tage 2009) to 
determine the opacity, the ratio of specific heats (7) and the mean 
molecular weight (/i) as a function of volume density and sound 
speed. 

This allows us to rapidly produce a suite of disc models span- 
ning a broad parameter space in M and rout • Figure [T] shows 
the required disc-to-star mass ratios for disc models with a given 
M and rout. It is important to note that quite high values of q 
are required to sustain self-gravitating discs - for typical T Tauri 
values (say M = lO'^M© j/r"\ rout = 100 au), the disc-to- 
star mass ratio required is g ~ 0.2 5, several orders of magni- 
tude higher than obser vations suggest dAndrews & Williamsl2007l : 
lOreaves & Ricell2010t) . If the local temperature exceeds 1000 K, 
then the ionisation fraction should be sufficiently high for MRI to 
operate, and enhance the local accretion rate. We model this by set- 
ting a — aMRi ~ 0.01, and re-calculate E such that the steady- 
state accretion rate is obtained. While this differs fr om the expected 
episodic nature of MR I activation in this scenario JArmitage et al.l 
uOOlMZhu et al.l2009f) , its draining effect on the resulting disc mass 



is important. For modest outer radii, the disc-to-star mass ratio in- 
creases to over 0.5 for all stellar masses. Even with MRI activa- 
tion providing drainage, the discs remain sufficiently cool that large 
amounts of matter can still be found in the inner regions (this is due 
to a decreasing steeply towards the central object). 

Models which satisfy the Tj criterion are classified as frag- 
menting models, and the local Jeans Mass is then calculated using 
equation ( I13t . We investigate three different stellar masses, gener- 
ating 250,000 models for each. The models span accretion rates 
from AI = 10^^ — lO^^A/0 yr~^, and disc outer radii from 
1.3 — 150 au (the inner disc radius is assumed to be 1 au for all 
models). 



5 RESULTS & DISCUSSION 

Figure [2] shows the expected Jeans Mass at fragmentation as a 
function of M and the disc outer radius rout- We present data 
for three different stellar masses: 0.5, 1 and 2Mq. The minimum 
Jeans mass across all stellar masses remains roughly constant at 
around 3 A/jup, as does the maximum Jeans Mass (~ 125 i\/jup, 
or O.IM0). These discs are defined to be quasi-steady with con- 
stant disc mass and local temperature - the minimum Jeans mass is 
therefore fixed by the heating/cooling balance in the disc (and by 
extension the local opacity). 

The minimum mass agrees well with iKratter et alj J201CI) 
(where they state a similar value for the minimum fragment mass at 
a temperature of 10 K and an accre tion rate of lO'^Mptyr"^), a s 
well as the results of lRafikovl |2005^ and'N ero & BiorkmanI ( l2009h . 
This demonstrates once more that objects formed by disc fragmen- 
tation will not occupy the Earth-mass regime (at least initially). If 
disc fragmentation is to produce such objects, the subsequent evo- 
lution of the fragments must include significant mass loss, as well 
as strong core sedimentation. While it has been shown that frag- 
ments of mass < IMjup are slow to pro duce substantial cores 
by sedimentation dUelled & Schuberlll2008l) , the typical fragments 
produced in these models may be more successful, especially with 
the dust-concentrating e ffects of the spi r al structure induced by 
gravitational inst ability ('Rice et al.' 2004; ' Clarke & Lodatol 120091 : 
lBoIev&Duriseni2010 : Nayakshin 2010b) . 

As has been found many times previously, discs with outer 
radii less than around 40 au will typically not form fragments 
(Rafikov 2005; Matzner & Levin 2005; Whitworth & Stamatello^ 
20061: lBolevetal.1 120061: iStamatellos & WhitworthI 12008 



Forgan et al.l |2009| ; Iclark3 |2009|) . While this is more a symp- 



tom of the imposed disc model (which in itself is not a particularly 
new implementation), we can use this result to satisfy ourselves 
that the model is indeed correctly distinguishing the fragmentation 
boundary. The dependence of the radial fragmentation boundary 
on stellar mass (i.e. lower mass stars have a lower boundary) is 
also an indication that the model is indeed performing as expected. 

The distribution of the contours shows a clear trend as the stel- 
lar mass is increased - the A/* — 0.5Mq star has contours more 
closely packed along the r-axis, and is able to form (for example) a 
20 Afjup fragment at 120 au (at the lowest M which permits frag- 
mentation at that radius). By contrast, the M, = 2Mq star can only 
form a 10 A/,jup fragment at the same location in M — r space. This 
is sensible, given that the disc around the A/* = 2A/q star has a 
much lower mass ratio for those model parameters. 

This can be seen in Figure[3] which shows the expected Jeans 
mass (in Jupiter Masses) at fragmentation, as a function of the disc- 
to-star mass ratio, q. We calculate models for which q is very much 
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Figure 1. 2D contour plots of the disc to star mass ratio, as a function of tlie steady state 'pseudo-viscous' accretion rate M, and the disc's outer radius rou 
for stellar masses of M* = O.SMq (top), M, = IMq (middle) and M* = 2Mq (bottom). 



higher than observed. Indeed, current mod elling indicates that q ^ 
1 disc s can only exist at very early times jMachida & Matsumotd 
|2010|). and they ev olve rapidly towards states where q < 0.5 
I Forgan et alj|201lh . The fact that typically q > 0.5 is required for 
fragmentation strains our assumption of local transport to breaking 
point. Further to this, requiring such large values of q further boosts 
the argume nt for rapid giant plan et formation while the system is 
still young JGreaves & Ricel2010l) . 

If we consider discs with g = 1 and rout ~ 50 au, then it 
becomes clear that the M« — 2Mq star can form the most mas- 
sive fragments (> 20M,jup), whereas the A/* — 0.5A1q star can 



only form fragments of mass ~ 7A/jup, with the M, — IMq star 
forming fragments of around lOA/jup. 

In the same fashion, the discs around the M* — 0.5Mq star 
require q ~ 1.25 to form a 10 A/jup fragment, while discs around 
the A/, = 2Mq require q ^ 0.75 (for rout = 50 au). We can 
identify a general trend - discs around higher mass stars can satisfy 
a maximum a fragmentation criterion (and hence form a minimum 
fragment mass of ~ 3Afjup) at lower q than discs around lower 
mass stars. 

What can be said about the expected statistics of these ob- 
jects? We present binned frequency distributions of the fragment 
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10-" 




Figure 2. 2D contour plots of the expected Jeans mass in spiral perturbations (in Jupiter masses), as a function of the steady state 'pseudo- viscous' accretion 
rate M, and the disc's outer radius Tout, for stellar masses of A/» = 0.5Mq (top), Af» = IMq (middle) and Af* = 2Mq (bottom). 



mass from the parameter survey in Figure|4](left panel). While this 
data will be in no way representative of the true distribution of frag- 
ment mass, it indicates that in a uniform disc population, certain 
fragment masses are considerably more popular than others (given 
the limits of the parameter space that we have investigated). In par- 
ticular, the peak of the distribution occurs between 10 and 20 i\f jup 
for all three stellar masses. Given that this is the regime in which 
the giant planet/brown dwarf boundary is located, this has inter- 
esting ramifications for planet formation theory, in particular for 

' " — ^5 

the canonical disc instability theory (Boss 1997), and th e related, 

more recent tidal downsizing theory (Nayakshin 2010a ; Bolev et al.l 



uOlOl) . Indeed, if this distribution is convolved with a more realis- 
tic distribution of disc parameters, this should provide an initial 
population of fragment masses. The onus is then on theory to be 
able to explain the evolution of this population into the current 
exoplanet population. Interestingly, the peak at 10-20 Mjup coin- 
cid es with the upper li mit for efficient core sedimentation derived 
bv iNavakshinI ( 1201 Obi) . If fragments are more commonly formed 
within the limits of efficient sedimentation, then it might be rea- 
sonable to assume that giant planets are cap able of hosting more 
substantial cores than previously assumed (cf lD'Angelo et al.l2010l 
and references within). 
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Figure 3. 2D contour plots of the expected Jeans mass in spiral perturbations (in Jupiter masses), as a function of the disc-to-star mass ratio q, and the disc's 
outer radius rout, for stellar masses of M* = 0.5Mq (top), M, = IMq (middle) and M« = 2Mq (bottom). 



Lower mass fragments occur more frequently around lower 
mass stars, which would suggest that if disc instability is to form 
planets, it will be required to do so at an early stage in protostel- 
lar evolution, as suggested by Greav es & Ricd ( 120101) . This is sup- 
ported by the frequency distribution of q (Figure ID right panel). 
Lower mass stars require in general a larger q to form fragments, 
again suggesting that giant planet formation by disc instability must 
occur early in the star's evolution. The minimum q for fragmen- 
tation is approxim ately 0.5, in agreement with the limit set by 
lKratteretai] ( l2010l) . 

We must now return to something of a thorny problem in mas- 
sive, self-gravitating discs - the a-parametrisation. We have demon- 



strated that fragmentation occurs when q > 0.5, and we have im- 
plicitly assumed a local model of angular momentum transport in 
our derivation of F./ and in the one-dimensional disc models. But 
such massive sel f-gravitating discs are not expected to be local 
JBalbus & Papaloizou.1999.) . The development of large amplitude, 
m = 2 spiral modes in the disc will produce non-local features 
in the gravitational potential, significantly affecting the stresses de- 
veloped at any location in the disc. Forganetal. (2011) demon- 
strate that non-local effects become important for isolated discs 
with q > 0.5, and also that Q and T become extremely vari- 
able with time. For discs with infall, this limit drops to g ^ 0.1 
dHarsono et alj|201 ih . The implication of this is that we should be 
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Figure 4. The frequency distribution of fragment mass (left) and disc-to-star mass ratio q (riglit) from disc models which resulted in fragmentation. Both 
distributions are binned - the bin widths are 0.2 Mj^p for the fragment mass, and 5 X 10^^ for q. 



cautious in applying the quantitative results from the disc models. 
To correct this, we require data on how A S/S depends on di sc 
variables when g > 0.5, in essence to repeat lCossins et alj J2009f) 's 
analysis at higher q. Equally, we must correct a similar deficiency in 
Fj by replacing the viscous heating term with something more ac- 
curate, perhaps also incorporating the effects of stellar irr adiation, 
an important s tabilising proces s in self-gravita ting discs jRafikovl 
l2007l : lcaiet"ai . 2008; Stamatel los & Whitworth 2008). If Q varies 
significantly, it may be prudent to investigate r./'s dependence on 
an extra timescale: 



= -^n 



(45) 



and recalculate Mj with the Q ~ 1 assumption relaxed. Indeed, 
its reliance on such few assumptions lends itself to use in numer- 
ical simulations where non-local transport does not dominate, and 
provides another means of ensuring that simulations satisfy the res- 
olution criteria of Bate & Burkert ( 1997), which explicitly requires 
the Jeans mass to be resolved by a sufficient number of SPH parti- 
cles (or equivalently that the Jeans length be resolved by a sufficient 
number of grid cells). 

While we have focused exclusively on self-gravitating proto- 
stellar discs in this paper, the Jeans mass expression permits itself 
to be taken to any size scale. For example, this work would apply 
to the fragmentation of discs around the supermassive black holes 
(SMBHs) that are expe cted to reside at the centre of most galaxies. 
iNavakshin et alj J200 7) use constant /3c SPH simulations to investi- 
gate star cluster formation (and subsequent star formation) around 
the Milky Way's own SMBH, Sgr A*. They find a different value 
for /3crit than is typically found for protostellar discs, and the local 
initial mass function (IMF) due to disc fragmentation is top-heavy 
in comparison to the Solar Neighbourhood. We find that (using 
their parameters of Afbii = 3.5 x IQ'^Mq, and Matsc = O.OlAfbh), 
that the Jeans mass is 

204 



Mj^ 



-Ma 



(46) 



(1 + i/VK) ' 

INavakshin et alj ( 120070 use /3c ~ 0.2 as their lowest value, giving 
a minimum Mj = 63Mq. This is still a round an order of mag- 
nitude larger than the masses obtained bv'Navaks hin et alj i2007t) 
and other works in the same field (e.g. Bonnell & R icell2008h . It is 
difficult to compare initial fragment masses with the resulting stel- 



lar masses grown through accretion, and it is also true that fragment 
destructi on and mass loss by fragment-fragment collision is more 
frequent JLevinl 20071) . Equally, this larger Jeans mass does suggest 
that the subsequent IMF will be top-heavy. An avenue for further 
work is to apply the expressions for Tj and Mj in the same fashion 
as we have done for protostellar discs to yield the initial distribution 
of fragment masses for a vari ety of AGN d isc parameters, and con- 
firm similar analyses by e.g. iLevinI J2007r) to obtain the expected 
mass distribution. 

We have said much about constant /3c simulations in this pa- 
per. While perhaps a little unrealistic, they have been extremely 
valuable to accretion disc researchers in getting a feel for what disc 
regimes are likely to produce fragments. Now that we have intro- 
duced a more generalised timescale F./, can we now consider a 
set of controlled experimental simulations where Fj is held con- 
stant? Mayer et al. (2004) perform a set of SPH simulations where 
Fe is controlled by allowing the mass of particles to increase by 
some prescription. While they use a more complex cooling func- 
tion than /3c, it is possible that the two techniques could be merged 
to produce a simulation where Tj is roughly constant at all radii, 
although incorporating a non-zero /3c is more challenging. Devel- 
oping such constant Tj simulations to gather empirical data is def- 
initely worth future investigation. 

We should finally acknowledge a subtlety regarding the sound 
speed. As with the surface density, we should really consider the 
sound speed inside the spiral wave, which will be enhanced relative 
to the background due to shocks. Shocks in self-gra vitating proto- 
stellar discs typically have Mach numbers M ~ 1 JCossins et al.l 
2009). The fractional perturbation amplitude in temperature AT/T 
in adiabatic shocks at 7 — 5/3forM — 1 is approximately 0.3125, 
implying that our expression for the Jeans Mass will underestimate 
the true value by roughly 30%. Preliminary work on Smoothed Par- 
ticle Hydrodynamics (SPH) simulations of disc formation from col- 
lapsing molecular clouds shows that the fragment masses are sys- 
tematically under-estimated by a similar factor (Forgan and Rice, 
in prep.). In any case, equation l ll3t is a minimum fragment mass 
estimate, and should be treated as such, especially when the growth 
of the fragments must also be considered. 

Indeed, an interesting avenue for further work is to follow the 
subsequent evolution of these fragments. Being able to accrete from 
what is still a substantial pool of gaseous material in the disc will 
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significantly influence tiieir final mass, and the available feedstock 
for any subsequent planet formation in the inner regions (say by 
core accretion). This will surely affect any potential bimodal dis- 
tribution of planetary parameters (Bolev 2009). As we have said 
previously, the ability of the fragments to grow as they potentially 
migrate inwards is of high importance both to canonical disc in- 
stability theory and to the more recent tidal downsizing hypothesis 
JNav akshin 2010b). Also, the future interaction of the fragments 
with the remaining disc will create strong tidal structures, provid- 
ing pressure maxima for dust grains to occupy and potentially en - 
hance their growth rate dRice et al.|[2004l : Iciarke & Lodatoll2009h . 
The expressions for Mj and Vj will be of use for further work in 
this two-phase regime. 



6 CONCLUSIONS 

We have derived the local Jeans mass inside a spiral perturbation, 
and shown that it can be calculable using azimuthally averaged 
disc variables. This allows the expected mass of objects formed 
from disc fragmentation to be calculated simply, both in theoreti- 
cal models and in numerical simulations, across a variety of size 
scales. We have taken this expression and derived a Jeans timescale 
r J, which describes the length of time required for the Jeans Mass 
to increase or decrease by its current value (normalised by the lo- 
cal angular velocity, in much the same way as the cooling time is 
often normalised). We show that the resulting expression for V j 
reduc es to the standard cooling time cr iterion as was previously de- 
fined ( lGanmiiell200ll ; lRice et al| 20051) . and fold s in related results 
on the disc's thermal histor y JClarke et afluOOTr) and envelope ac- 
cretion JKratter et akluOlOf) . We confirm that subjecting the disc to 
extra shock heating promotes fragmentation, and that rapidly ac- 
cretion encourages discs to fragment, provided that local angular 
momentum transport and cooling is efficient. 

We have investigated the functional form of the Jeans Mass by 
applying it to simple one-dimensional steady-state self-gravitating 
disc models. We find that the minimum fragment mass is around 
3 Jupiter Masses, in agreement to previous work on the su bject 
( lRafikov|[2003 : iKratter et ai]|2010l ; Ffero & Biorkmarj|2009h . We 
also find that in a uniform survey of disc parameters, the most fre- 
quent fragment masses lie between 10 - 20 A/jup. These results 
are relatively insensitive to stellar mass, although lower mass stars 
require a higher disc-to-star mass ratio to produce fragments than 
their higher-mass counterparts. Conversely, the mass of fragments 
found around lower mass stars are higher than that for higher mass 
stars (for the same outer radius and steady-state accretion rate). We 
do however note that the derivation of the Jeans mass relies on as- 
suming local angular momentum transport, which is often not the 
case for discs of high disc-to-star mass ratio. The errors introduced 
by this assumption require further empirical data about non-local 
angular momentum transport to quantify accurately. 

The formalisms introduced in this paper are amenable to a va- 
riety of studies in self-gravitating discs of all scales, allowing the 
intial fragment mass distribution to be identified given the known 
disc population, and to provide more general criteria for disc frag- 
mentation itself. 
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